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SUMMARY 

A  source-distribution  method  is  applied  to  derive  a  solution  for 
the  time-dependent  surface  velocity  potential  of  thin  finite  wings  at 
supersonic  speeds.  The  solution  is  illustrated  by  evaluating  the 
upwash  over  the  tip  of  an  arbitrary-plan-boundary  wing  having  a 
supersonic  and  subsonic  leading  edge.  The  upwash  then  gives  the 
effective  sources  of  the  flow  field  lying  between  the  wing  plan 
boundary  and  the  foremost  Mach  waves,  which  are  applied  in  the 
derivation  of  the  wing-surface  velocity  potential.  A  simple  example, 
the  load  distribution  on  a  wing  whose  effective  angle  of  attack  is 
changing  linearly  with  time,  is  included  to  illustrate  the  appli¬ 
cations  of  the  derived  expressions. 


INTRODUCTION 

The  analysis  of  the  aerodynamic  effects  in  the  vicinity  of  thin 
wings  at  supersonic  speeds  can  be  simplified  to  obtain  useful  results 
by  means  of  the  linearized  theory.  Steady-state  or  time- independent 
solutions  have  been  obtained  for  enough  cases  (for  example,  refer¬ 
ences  1  to  27 )  that  the  essential  features  of  the  load  distributions 
on  various  parts  of  the  wing  can  be  determined  by  graphical  or 
analytical  methods.  The  time-dependent  load  distributions  are  more 
difficult  to  obtain.  Such  problems  include  the  transient  effects 
of  gusts,  changes  in  angle  of  attack,  skin  vibration,  and  flutter. 

A  number  of  investigators  have  studied  two-dimensional  time- 

dependent  flows  over  thin  wings.  These  flows  are  generally  included 
as  special  cases  of  the  theory  of  reference  27.  The  method  of 
reference  27  is  similar  to  the  steady-state  solution  of  reference  1 
and  includes  three-dimensional  or  finite  wing  solutions  in  cases 
where  the  aerodynamic  effects  of  the  bottom  and  top  wing  surfaces  are 
independent.  No  solutions  are  known  to  have  been  published  for 
cases  involving  interaction  between  the  flow  over  the  bottom  and 
top  wing  surfaces. 
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The  concepts  of  references  9  and  18  that  lead  to  solutions, 
which  include  the  effects  of  interaction  between  the  bottom  and 
top  wing  surfaces  in  the  steady  state,  can  also  be  applied  to 
obtain  time-dependent  solutions.  The  derivation  of  the  surface 
velocity  potential  in  regions  influenced  by  subsonic  leading 
edges  was  completed  during  January  1948  at  the  NACA  Cleveland 
laboratory  and  is  presented  herein. 


ANALYSIS 

In  order  to  unify  the  discussion,  parts  of  the  fundamental 
treatment  presented  in  reference  27  are  repeated.  The  analysis 
includes  the  derivations  of  (l)  the  time-dependent,  linearized 
partial-differential  equation  for  the  perturbation-velocity 
potential  of  an  ideal  fluid,  (2)  the  fundamental  solution  that 
will  satisfy  the  boundary  conditions  on  the  wing,  (3)  the  upwash 
between  the  wing  boundary  and  the  foremost  Mach  line,  and  (4)  the 
velocity  potential  on  the  surface  of  the  wing. 

Differential  equation.  -  The  linearized  Euler's  equations 
for  a  compressible  fluid  may  be  written 


+U^P.  _ 

dxdt  bx2  pQ  bx 

d2<P  1  dp 

dydt  <5x5y  *  c5y 

1  dp 

dzdt  dxdz  Pq  bz 

where 

cp  perturbation-velocity  potential 


t  time 

U  free-stream  velocity 

p^  free-stream  density 


(1) 
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* 


p  static  pressure 

/ 

x,y,z  Cartesian  coordinates  (free  stream  parallel  to  x  axis) 

(A  complete  list  of  symbols  is  included  in  appendix  A.) 

When  equations (l)  are  multiplied  by  dx,  dy,  and  dz  respec¬ 
tively,  added,  and  integrated,  the  result  is 

f: +  °r  *  *  -  «<*>  <2> 

at  ax  pQ 

where  g(t)  is  an  integration  constant  at  any  given  time. 

The  linearized  continuity  equation  can  be  written 


at  ax  0  Vax2  ay2  az2/ 

or,  because  the  velocity  of  sound  c  is 


°2-£~o2 


equation  (3)  becomes 


(jL  +  ugA  ,e  afj  a^p  afjp 

v5*  P0c02  ay2  az2  0 


(3) 


(3a) 


Substitution  of  p/p0  from  equation  (2)  in  equation  (3)  yields 


(1  .  „2)  is  +  is  +  £2  .  i  +  22  +  i  (4) 

dx2  dj2  Si2  o2  at2  o2  Sx&t  c2  at 

where  M  is  the  Mach  number  of  the  free  stream  and  the  zero  sub¬ 
script  has  been  dropped  from  c. 

Equation  (4)  is  the  required  linearized  partial-differential 
equation  for  the  velocity  potential.  If  cp  is  independent  of  time, 
the  Prandtl-Glauert  equation  immediately  results. 
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The  function  g(t)  depends  on  the  condition  of  the  flow  ahead 
of  the  body.  If  the  flow  Is  uniform  and  undisturbed,  g  will  be 
constant  (from  equation  (2))  and  equal  to  Pq/Pq*  Tk®  function  g 

will  be  assumed  constant  in  the  rest  of  the  analysis. 

A  change  In  variables  will  convert  equation  (4)  to  a  standard 
form  of  the  wave  equation.  The  transformation  equations  are 


Basio  solutions  of  equation  (6)  corresponding  to  spherical  waves 
are 


where  r’  ■  yi’2  +  y'2  +  z'2,  02  =  M2  -  1,  and  f  is  an  arbitrary 

function. 


The  basic  solution  for  the  supersonic  case  Is  obtained  as  the 

sum  of  equations  (7).  (See  reference  27.)  If  this  solution  Is 
transformed  to  a  general  point  in  the  x,y,t  space,  the  basic 
solution  of  equation  (4)  assumes  the  form 
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where  w  Is  the  local  z  component  of  the  perturbation  velocity. 
Equation  (11)  defines  the  strength  of  the  local  source  in  terms  of 
the  z  component  of  the  perturbation  velocity.  The  fact  that  equa¬ 
tion  (ll)  is  independent  of  the  value  of  t ^  is  in  agreement  with 
the  statement  that  only  those  subelements  near  the  point  (x,y,0) 
contribute  to  the  vertical  velocity  Scq/Sz  on  the  wing  surface. 
Physically,  this  evaluation  infers  that  the  flow  is  tangent  to  the 
wing  surface  at  each  local  point. 

The  time-delay  terms  of  equation  (8)  can  be  denoted  by  Ta  and 
T  ^  where 


The  velocity  potential  at  any  point  (x,y )  is  obtained  by  integrating 
the  sources  in  the  z  =  0  plane  over  the  area  S  included  in  the 
forward  Mach  cone.  By  use  of  equation  (11),  the  velocity  potential 
becomes 


°  [w(M*t-Ta>+w(  ^n,t-Tb)]  d^m 

s  ^(x-D2  -  P2(y-!i)2  -  P2z2 

Equation  (13)  was  derived  by  Garrick  and  Rubinow  in  reference  27  and 
gives  the  velocity  potential  at  any  point  in  space  in  terras  of  the 
z  component  of  the  perturbation  velocity  in  the  z  =  0  plane.  The 
equation  reduces  to  the  steady-state  solution  of  Puckett  (reference  l) 
when  w  is  independent  of  time. 

/ 

A  physical  interpretation  may  he  given  for  the  time  delays  T 

a 

and  of  equations  (12)  and  (13).  If  a  disturbance  is  generated 

at  point  (|,T))  at  time  t  =  0,  the  wave  front  from  that  disturbance 
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will  travel  outward  as  a  spherical  wave  about  a  center  that  moves  with 
the  free- stream  velocity.  (The  trace  of  the  waves  on  the  z  ■  0  plane 
is  illustrated  in  fig.  2.)  The  wave  front  will  enter  and  emerge  from 
the  point  (x,y,z)  at  two  later  times  and  T^.  The  equation  of 

the  spherical  wave  path  that  passes  through  the  point  (x,y,z)  is 


(x-l-TJT)2  +  (y-q)2  +  Z2  m  c2!2  (14) 

The  solution  for  T  is 

T  .  <x4)m+  ijxAf.  -  p2 (y-ri)2  -  P2*? 

P2c  p2c 

which  gives  equations  (12).  (At  a  given  point  (x,y,z),  the  strengths 
of  the  same  wave  at  the  two  times  Ta  and  are  equal  despite  the 

change  in  the  radius  of  the  wave  front.)  At  a  given  time  t,  only  the 
wave  fronts  that  are  entering  and  emerging  from  the  point  (x,y,z)  con¬ 
tribute  to  the  velocity  potential.  These  two  waves  originated  at  the 
point  (|,t],0)  at  times  (t-Ta)  and  (t-Tb). 

The  remainder  of  the  analysis  is  primarily  concerned  with  the 
aerodynamics  in  the  plane  of  the  wing  so  that  z  may  be  set  equal  to 
zero.  Equations  (12)  and  (13)  may  also  be  conveniently  expressed  in 
an  oblique  coordinate  system  whose  axes  lie  parallel  to  the  Mach  lines. 
(See  fig.  3.)  The  transformation  equations  are 


v  =2L  (Upr,) 

2p 

e  - 1  <t«) 

n  =  i  (v-u) 

“» -  §  (l“P7) 

vw  “  2p 

1  ■  |  <V"w> 

7  “  M  (vVUw) 

(15) 

Inasmuch  as  the  elemental  area  in  the  (u,v)  coordinate  system  is 

dudv,  equations  (13)  and  (12),  in  the  case  of  z  *  0,  become 
Mz 
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q> 


2M« 


'S 


[w(u,v,t-Ta)  +  w(u,v,t-Tb)  dudv 

V(uw“u) 


(16) 


M(yw-v+uw-u)  +  2  YT»V-u)(^-y) 
a  Mpc 


T 

b 


M(vv-v+\~u)  -  2V(Vu)(yw-y) 

Mpc 


(17) 


Equation  (16)  gives  the  velocity  potential  in  the  t  «  0  plane 
in  terms  of  the  perturbation-velocity  component  v  normal  to  the 
plane.  If  only  supersonic  leading  edges  are  included  in  the  forward 
Mach  cone  from  (x,y),  w  may  be  evaluated  in  terms  of  the  effective 
wing  slopes  O  measured  in  t)  *  constant  planes  by  the  relation 

w  -  U a  (18) 

If  a  subsonic  leading  (or  trailing)  edge  is  also  included  in  the 
forward  Mach  cone  from  (x,y),  the  slopes  of  the  stream-lines  X 
associated  with  the  upwash  between  the  wing  boundary  and  the  foremost 
Mach  wave  must  be  evaluated  and  included  in  the  calculation  (equa¬ 
tion  (16))  for  the  velocity  potential. 

Upwash  between  wing  boundary  and  formost  Mach  line.  -  The  slopes 
of  the  streamlines  in  the  region  Sj>  of  figure  4  could  conceivably 
be  generated  by  a  thin  wing  or  diaphragm,  as  employed  in  reference  9. 

The  flows  above  and  below  the  z  *  0  plane  may  then  be  independently 
treated.  The  velocity  potential  on  the  top  surface  of  the  z  *  0 
plane  is  given  by  equations  (16)  and  (18)  as 
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CD _ U_  /  /  EtK^-V  +aT(u>v>t'V 

t  ■  as;  I  I  - 

yCup-uHYp-v) 


dudv 


_U 

2Mir 


[x(u,v,t-Ta)  +  7<u,Y,t-Tb)^|  dudv 

s  y(uD-u)(vD-v) 


(19) 


where  Up  and  Vj,  are  the  coordinates  of  the  point  at  which  cp  la 
evaluated,  0^,  represents  the  slopes  of  the  streamlines  on  the  top 
wing  surface,  and  A  represents  the  slopes  of  the  streamlines  in  the 
field  Sp  from  the  point  of  view  of  the  top  wing  surface.  Similarly, 
the  potential  on  the  "bottom  surface  of  the  z  =  0  plane  is 


0B(u,v,t-Ta)  +  O^UjVj^T^  dudv 
/^(ud-u)(vd-y) 


u,v,t-Ta)  +  A(u,v,t-Tb)J  dudv  ^  ^ 

/|/(urru)(vI)-v) 

The  pressure  at  a  given  point  of  the  field  SD  (fig.  4)  In  the 
plane  of  the  wing  can  be  calculated  by  substituting  either  cpT  or 
cpg  Into  equation  (2).  The  two  computations  of  the  pressure  can 
then  be  equated: 
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Equation  (21)  has  the  solution 

CPT  -qpfi  +  2E(x-Ut,y)  (21a) 

where  E  is  an  Integration  function.  Substitution  of  equations  (19) 
and  (20)  In  (21a)  yields 


,-JL  f  f  <V°r>..i,dalT  tH 

Sv 

(21b) 


(In  equation  (21b)  and  In  some  of  the  equations  to  follow ,  the 
notation  Xft  means  X  (u,v,t-Ta)  and  X  a  ^  means 

|~X(u,v,t-Ta)  +  X(u,Y,t-Tj))J.  A  similar  notation  is  applied  to 

0B>  Ot,  and  (o^-Op).)  If  aquation  (21b)  Is  substituted  Into  (19), 

the  result  Is 


(W..bdudT 

2  ^(ub-«)(td-t) 


+  H 


(22) 


Equation  (22)  represents  the  velocity  potential  In  the  plane  of  the 
wing  for  the  region  Sp.  Similarly,  formulation  of  would  change 

only  the  sign  affixed  to  H.  The  function  2H  of  equation  (21a) 
represents  the  difference  In  potential  across  the  z  «  0  plane, 
corresponding  to  the  strength  of  vorticlty  in  the  wake  of  the  wing 
(reference  18).  For  flat-plate  wings,  Og  +  0ip  »  0  and  H  Is 
Just  the  potential  on  the  top  surface  of  the  vortex  sheet. 

The  foremost  Mach  wave  (fig.  4  or  5)  originating  on  the  leading 
edge  generally  represents  a  line  of  Infinitesimal  disturbance  along 
which  E(x-Ut,y)  can  be  set  equal  to  taro  at  all  times.  The  func¬ 
tion  E  remains  zero  along  y  «  constant  lines  for  values  of  x 
not  intercepted  by  the  wing  or  a  material  body  (region  Sp  ^  of 
fig.  5).  The  region  Sp  2  of  figure  5  generally  contains  a  vortex 
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sheet  lying  in  the  plane  of  the  wing  and  H  is  not  zero.  The  func¬ 
tion  H(y),  established  along  the  wing  trailing  edge  at  some  time  t, 
remains  unaltered  for  later  times  along  a  curve  that  sweeps  down¬ 
stream  with  the  free-stream  velocity  and  has  the  form  of  the  wing 
trailing  edge.  The  rest  of  the  discussion  is  concerned  with  only  the 
effects  of  leading  edges,  that  is,  those  cases  for  which  H  is  zero. 

The  origin  of  coordinates  for  the  wing  shown  on  figure  4  is 
placed  at  the  junction  of  the  supersonic  and  subsonic  leading  edge. 

The  supersonic  and  subsonic  leading  edges  are  respectively  defined  by 
the  equations 


*5? 

rH 

> 

II 

> 

or 

u  =  u^(v) 

v  =  v2(u) 

or 

u  =  u^v) 

(23) 

Equation  (21b)  with  H 

set  equal  to  zero, 

then  becomes 

(21c) 


In  the  steady-state  solution,  X  and  (cJg-Oip)  are  independent 
of  Up  (references  9  and  ll),  so  that  the  integrations  with  respect 
to  v  may  be  equated  to  give 


The  reduction  of  equation  (21c)  to  give  equation  (24)  for  the 
time- independent  case  shows  that  only  those  wing  slopes  along 
u  =  up  contribute  to  the  upwash  for  points  on  this  Mach  line.  In 

the  time-dependent  case,  both  \  and  (Og-Cj)  contain  functions 
of  ud,  and  the  validity  of  the  reduction  may  be  questioned.  If 
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in  the  time-dependent  case  only  those  wing  slopes  along  u  =  uD 

were  to  contribute  to  the  upwash  for  points  on  this  Mach  line, 

X  and  (Op-dp)  would  become  independent  of  up  and  the  reduction 
of  equation  (21c)  to  equation  (24)  would  again  be  Justified.  An 
argument  to  substantiate  this  choice  is  presented  in  appendix  B. 

If  u  is  replaced  by  uD  in  equation  (24),  the  two  time 
delays  become  equal  and  are  linear  with  respect  to  y: 


Each  ininitesimal  wing  element  then  produces  at  time  t  and  place  v^ 
an  increment  in  upwash  corresponding  to  the  steady-state  effect  of 
wing  elements  whose  slopes  are  the  same  as  the  time-dependent  elements 
evaluated  at  time  t-Ta.  The  solution  of  equation  (24)  is  the  sum 

of-  a  series  of  infinitesimal  steady-state  solutions,  each  of  which 
satisfies  Abel's  equation  (reference  11)  and  each  of  which  requires 
integration  only  over  the  wing  slopes  along  u  =  uD. 

In  order  to  illustrate  the  argument,  the  increment  in  X  at 
point  (up,  Vp)  (fig.  6)  due  to  a  steady-state  wing  element  of 

length  dv  at  point  (uD,v)  may  be  obtained  (appendix  B)  from 
Abel’s  equation  in  the  manner  of  reference  11  as 


6.X  (°B’0T^  iv2~v 

dv  =  _ 

2jt(vD-v)/^  vD-v2 


(25) 


Equation  (25)  may  be  applied  in  either  the  time-dependent  or  the 
time- independent  cases.  For  the  time- independent  cases,  the  wing 
slopes  are  evaluated  at  uj),v.  For  the  time-dependent  solutions, 
the  wing  slopes  are  evaluated  at  (u^v)  at  time  t  - 

Integration  of  equation  (25)  across  the  wing  gives 
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A(u,vr.,t)  = 


2« 


VvD'v2 


0v< 


v^-v 


vD-y 


(26) 


where  0-g  T  means  ((7^-0^) .  That  equation  (26)  is  a  solution  of  equa¬ 
tion  (24)  is  shown  in  appendix  B. 

Evaluation  of  velocity  potential  on  wing  surface.  -  Equations  (19) 
and  (26)  now  allow  the  calculations  of  the  velocity  potential  on  the 
wing  surface.  With  reference  to  figure  7,  the  velocity  potential  at 
the  point  (uw,vw)  is 


<P. 


7 T  .  dudv 
a.b _ 

Y(uw-u)(vw-v) 


U 

2Mrt 


n  n 


a.b 


dudv 


v,l+2 


U 


SD 


][(%- u)(vv-v) 
(27) 


For  the  wing  of  figure  7,  the  integration  of  X  over  the 


area  is 


* 


I  -  - 


U 

2M jt 


U 


X  dudv 
a.b _ 

/\/(uw-u)(vw-v) 


U 

2Mn 


U 


nu2<vw> 


pv 


du 


V(u„-u) 

'  u 


X  .  dv 
a.b  D 


fvvD 


(28) 


If  X  from  equation  (26)  is  substituted  into  equation  (28),  there 
results 
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Equation  (28b)  represents  the  contribution  to  the  velocity  potential  of  the  upwash  over 
the  wing  tip  and  can  be  used  to  rewrite  parts  of  equation  (27).  If  the  velocity  potential 
from  equation  (27)  is  written  in  full  for  the  wing  of  figure  7 ,  the  following  results  is 
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(continued  on  the  following  page) 
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Although  equation  (27a)  looks  rather  formidable,  there  are 
actually  only  two  types  of  term.  Simple  variations  are  obtained 
by  changing  signs  from  minus  to  plus,  by  replacing  bottom-surface 
wing  slopes  by  top-surface  wing  slopes,  or  by  altering  the  limits 
of  Integration.  Simplifications  also  occur  In  specific  examples, 
although  the  Indicated  Integrations  to  obtain  explicit  solutions 
are  generally  difficult  to  perform. 

If  the  wing  slopes  In  the  vicinity  of  the  wing  tip  vary  con¬ 
tinuously  with  either  time  or  position,  equation  (27a)  can  be 
applied  in  Its  present  form  without  altering  the  limits  of  Integra¬ 
tion.  The  equations  for  the  wing  slope  would  be  Inserted  Into 
equation  (27a)  and  the  indicated  integrations  conducted.  This  type 
of  calculation  has  been  Illustrated  for  periodic  oscillations  of 
flat-plate  wings  in  references  27  and  28  when  the  bottom-  and  top- 
wing  slopes  are  independent.  A  less  complicated  example.  In  whioh 
the  angle  of  attack  of  all  wing  elements  of  a  flat  plate  varies 
linearly  with  time,  illustrates  this  type  of  calculation. 

If  the  angle  of  attack  of  all  wing  elements  is  changed  at  a 
uniform  rate  m  for  all  times,  the  wing  slopes  may  be  expressed 
as 

0g  a  01  +  mb  m  —  O rp  (29) 

where  a  Is  the  angle  of  attack  of  the  wing  at  time  t  »  0. 

These  effective  wing  slopes  may  correspond  to  a  constant  accelera¬ 
tion  of  mU  in  the  z  direction.  Substitution  of  equation  (29) 
into  equation  (27a)  gives 


(30) 
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Integration  with  respect  to  v  produces 

PS 


q>„«22 

T  Mrt 


(a+nrt) 


u 


«2 


Vs-Ti<a> 

Vv* 


du  -  -5L 

3pc 


S 


S 


(▼y-Tl(tt)) 


du 


m 

pc 


PS 


s 


(30a) 


If  the  equation  for  the  leading  edge  of  the  wing  is  v  =  -  k^u, 
equation  (30a)  may  he  integrated  to  obtain 


9, 


_  2U 
T  Mr 


! 


a  +  nrt  -  JjL  (^  +1)^  +£(  ±;-§W 


2L/  i_-5 

4pcl  kx  3r* 


+  -JLj  i-  !iL 
2pc\  3/S 


V  <S“S^S+klS) 


(a  +  mt)(k]Uw  +vw)  - 


t 

(k1  +  l)a(k1uw  +vw) 


4pcki 


tan 


-1  iP'M) 


V  S  +klS 

(31) 

The  pressure  coefficient  nay  be  derived  from  equation  (2)  as 


■  .222. 

'P  x  2 

I°on 


S  S  +  4 


(32) 


Substitution  of  equation  (31)  into  equation  (32)  results  in  the  pres¬ 
sure  coefficient  on  a  family  of  wings  in  the  region  influenced  by 
the  subsonic  leading  edge. 


20  NACA  TN  No.  1699 

Page  20:  Equation  (33)  should  read  as  follows: 


|  WJ.CN 


(33) 

where  Ug(v)  Is  evaluated  at  v  =  vy.  (Although  equation  (33)  was 

derived  for  the  case  of  constant  acceleration  mU  in  the  z  direc¬ 
tion,  the  solution  may  he  combined  with  other  equations  to  evaluate 
the  pressure  distributions  for  a  variety  of  wing  motions.  For 
example,  a  wing  rotating  with  constant  rate  m  about  an  axis  of 
pitch  fixed  with  respect  to  the  wing  and  lying  in  the  z  *  0  plane 
would  have  the  pressure  distribution  given  by  equation  (33)  superposed 
on  the  pressure  distribution  associated  with  uniform  rates  of  pitch 
described  in  reference  21. ) 

The  steady-state  solution  included  in  reference  18  is  obtained 
from  equation  (33)  by  setting  m  =  0.  The  solution  for  the  infinite 
swept  wing  is  obtained  by  setting  vw  «  Ug(vv)  =  0  and  computing 
Cp.  The  values  of  Cp  obtained  along  the  line  vw  =  0  are  then 
constant  along  lines  parallel  to  the  leading  edge.  The  load  dis¬ 
tributions  of  a  family  of  wing  plan  boundaries  may  be  evaluated  by 
choosing  the  desired  equation  u  =  u2(v)  for  the  wing-tip  plan 
boundary. 

The  general  equation  (27a)  also  includes  solutions  for  which 
a  finite  number  of  discontinuities  can  exist  in  the  wing  slopes  with 
respect  to  either  position  or  time.  The  procedure  is  the  same  as 
in  the  case  of  continuous  wing  slopes,  except  that  the  fields 
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of  integration  will  be  subdivided  in  accordance  with  the  requirement 
of  the  discontinuity.  For  example,  in  the  problem  of  determining 
gust  loads,  the  effective  angle  of  attack  of  the  wing  will  change 
discontinuously  along  a  line  parallel  to  the  y  axis  that  moves 
downstream  with  the  free- stream  velocity.  If  the  initial  wing  slope 
is  taken  as  zero  and  the  increment  in  wing  slope  due  to  the  gust  is 
a,  the  slopes  of  the  wing  at  time  t  -  T  (as  time  appears  in  equa¬ 
tion  (27a))  are  either  zero  or  a  according  to  the  inequalities 

I  >  U(t  -t),  a£  =  -oT  =  o 

i  <U(t-T),  Oj  =  -Oij,  *  a 

The  two-dimensional  solution  of  this  problem  is  presented  in  ref¬ 
erence  29.  The  three-dimensional  solution  of  this  and  other  transient 
problems  within  the  scope  of  equation  (27a)  requires  further  research. 


Flight  Propulsion  Research  Laboratory, 

National  Advisory  Committee  for  Aeronautics, 
Cleveland,  Ohio  ,  May  18,  1948. 
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APPENDIX  A 
SYMBOLS 

The  following  symbols  are  used  in  this  report, 
pressure  coefficient, 

velocity  of  sound 
function 

integration  function  of  time  (herein  considered  as 
constant ) 

integration  function  of  i-Ut  and  y 
integral 

constant  greater  than  zero 

free- stream  Mach  number 

rate  of  increase  of  wing  angle  of  attack 

static  pressure 

^(x-  |)2  -  p2(y  -  n)2  -  P2z2 

+  y'2  +  *’2 

plan- form  area 
manipulation  variable 
time 

(l-M^  t  +2! 

c 

free-stream  velocity,  taken  parallel  to  x  axis 

oblique  coordinates  whose  azea  lie  parallel  to  Mach 
lines  in  z  =  0  plane 


P-Pp 
^  pu2 
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v  *  component  of  perturbation  velocity,  (taken  as 

positive  in  direction  of  an  outwardly  drawn  normal 
from  plane  of  wing) 

x.y,z  Cartesian  coordinates 


x’jy'jZ’  transformed  Cartesian  coordinates _ 

(x'  =  x,  y'  =  </l-M2y,  z'  ■  </l-M2z) 

|,q,£  Cartesian  coordinates  In  x,y,z  directions,  respectively 


a 

3 


angle  of  attack 


cotagent  of  Mach  angle. 


\  effective  slopes  of  streamlines  (measured  in  t)  *  constant 

planes)  in  z  =  0  plane  between  wing  boundary  and 

foremost  Mach  line,  X  =  ^ 


p  density 

a  effective  wing-section  slopes  measured  in  t)  =  constant 

planes  (  6  =  w/u) 

T  time  delay 

i 

cp  perturbation-velocity  potential 

Subscripts: 

0  free  stream 

1,2,3  .  .  .  numbered  areas  or  wing-plan-boundary  equations 
a,b  time  delays  Tft  and  T-^ 

B  bottom  (of  wing) 

D  upwash  field 

T  top  (of  wing) 


w 


wing 
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Examples: 


”1 

u2 


S 


w(l  +  2) 


du» 


dv. 


w 


Aa,b 

(cfe  -aT)a>b 


curve  v  =  v-^(u)  along  supersonic  leading  edge 
curve  u  =  Ug(v)  along  subsonic  leading  edge 
wing  areas  1  plus  2 

derivative  of  curve  U2(vw)  with  respect  to  v^ 

slope  at  time  t-  T&  plus  slope  at  time  t  -Tb 

difference  between  bottom  and  top  wing  slopes  at 
time  t  -Ta  plus  this  difference  at  time  t-Tb 


aB,T  aB  ”  °T 

u^,v^  oblique  coordinates  of  point  x,y  on  wing 
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APPENDIX  B 

EVALUATION  OF  UPWASH 

Origin  of  upwaah.  -  The  fundamental  source  solution  (equa¬ 
tion  (8);  for  the  three-dimensional  velocity  potential  of  thin 
wings  moving  at  supersonic  speeds  can  he  written  as 

cp  =  (Bl) 

If  this  equation  is  partially  differentiated  with  respect  to  z, 
the  result  is 

§2  _  fir  (b2) 

dz  dr  dz  r  dr 

If  z  is  allowed  to  approach  zero,  the  quantity  dcp/dz  will 
approach  zero  unless  i  approaches  infinity  (that  is,  r 

r  or 

approaches  zero).  The  fundamental  solution  (Bl)  can  then  he  integrated 
over  a  surface  of  sources  in  the  z  *  0  plane  to  obtain  an  extended 
velocity  potential.  Because  of  the  form  of  equation  (B2),  the  z  com¬ 
ponent  of  the  perturbation  velocity  generated  by  this  potential  will 
arise  only  from  those  subelements  near  r  =  0. 

On  the  wing,  the  flow  must  be  tangent  at  each  point  to  a  defined 
surface.  The  quantity  dq^dz  is  therefore  determined  at  each  local 
point  by  the  wing.  Under  such  circumstances  and  because  of  the  wing 
restraint,  only  the  subelements  near  r  «  0  obtained  by  assuming 
(x-|)«(y-Tj)=rf(z-t)  ^0  contribute  to  dcp/dz .  No  such  restraint 
exists  for  the  upwash  field.  If  the  same  condition  is  imposed  to 
evaluate  dq/dz  in  the  upwash  field  part  of  the  z  «  0  plane,  the 
trivial  answer  results, 

dz  dz 

A  fundamental  distinction  exists  between  the  evaluation  of 
dq/dz  on  the  surface  of  the  wing  and  in  the  upwash  field.  The 
wing  surface  is  a  restraint  and  can  generate  primary  impulses  that 
lead  to  the  velocity  potential  and  hence  to  the  flow  distributions. 

The  upwash  field  is  unrestrained,  however,  and  can  transfer  from 
the  bottom  to  the  top  wing  surfaces  only  those  impulses  that  have 
already  been  generated  by  the  wing. 
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The  upwash  must  arise,  however,  from  those  subelements  near 
r  =  0.  Because  the  unrestrained  subelementa  in  the  vicinity  of  the 
point  (x-f)  =  (y-T])=z  =  0  cannot  generate  upwash,  dop/dz 
(for  z  =  0 )  must  be  generated  by  those  subelements  in  the  vicinity 

of  the  curve  (x-|)2  -  02(y-f))2  =  02z2  =  0.  This  condition  implies 
that  only  those  wing  slopes  along  u  =  Up  generate  the  upwash  for 
points  on  this  Mach  line.  The  defining  equation  for  A  should  then 
be  restricted  by  setting  u  =  Up,  for  either  time- independent  or 
time-dependent  flows.  This  result  was  proved  directly  from  the 
integral  equation  (21c)  for  time- independent  cases  (reference  11). 

Differential  equation  for  A,  -  The  increment  dA  at  point 
Ujj,Vj)  due  to  the  wing  elements  dv  (fig.  6)  of  wing  slopes 

— r —  may  be  determined  from  equation  (6)  of  reference  11  as 


dA  s 


1  _d_ 

s  ds 


U 


ns 


dVi 


nv  +  dv 


(^  -  aT)dv 

2\N  • T 


(B5) 


where  s  is  the  manipulation  variable  and  is  evaluated  at  s  *  v^. 
The  integration  of  dv  over  an  infinitesimal  distance  results 
in  the  removal  of  one  integration  sign.  Also,  (<Xg  -  0<p)  may  be 
considered  as  constant.  Equation  (B3)  then  becomes 


dA 

dv 


gB~  gT 
2s 


b_ 

ds 


ns 


v2 


dv. 


-  vd)(td  " 


(OB  “  gT>  Vv2  -  Y 
2s  (s  -  v)^s  -  v2 


(B3a) 


Replacement  of  the  manipulation  variable  s  by  v^  gives  equa¬ 
tion  (25) 


dA 

dv 


(°B  " 

2it 


(td  -  v)  »/vD  -  ,2 


(B3b) 
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equation  (24)  when  u  *  uj,.  -  If  the  integration  variable  v  is 

changed  to  s  to  prevent  confusion,  equation  (26)  assumes  the 
following  form: 


The  first  member  of  equation  (24)  then  becomes 


r» 

_ dv _  _ _ « _ 

Jv  <v-s)  V(V_V2)(VD'V)  V^T2“8HtD-8) 
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Therefor© 


Equation  (B7 )  is  equation  (24)  for  u  *  u^. 
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Figure  1.  -  Field  of  integration  for  evaluating  velocity  potential 

(equation  (10))  of  a  source. 
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Figure  2.  -  Relation  between  time  delays  Ta  and  and 

position  of  wave  front. 
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Figure  3.  -  Comparison  of  Cartesian  and  oblique  coordinate  systems 


Figure  4.  -  Fields  of  integration  for  evaluating  upwash  between 
wing  boundary  and  foremost  Mach  wave* 
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Figure  7»  -  Integration  limits  for  equations  (27)  and  (28) 


